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Abstract 

We investigate the mass spectrum of l"*"^ exotic mesons with quark 
content (cqcq), using molecular and diquark-antidiquark operators, in 
quenched lattice QCD with exact chiral symmetry. For the molecular 
operator {(q7jc)(c75q) — (c7jq)(q75c)} and the diquark-antidiquark op- 
erator {(q-^C7ic)(qC75C-^) — (qC7j^c-^)(q^C75c)}, both detect a reso- 
nance with mass around 3890 it 30 MeV in the limit — > m^, which 
is naturally identified with X(3872). Further, heavier exotic meson res- 
onance with J^*-" = l"*""*" is also detected, with quark content (cscs) 
around 4100 ± 50 MeV. 



PACS numbers: ll.15.Ha, 11.30.Rd, 12.38.Gc, 14.40.Lb, 14.40.Gx 
Keywords: Lattice QCD, Exact Chiral Symmetry, Exotic mesons. 
Charm Mesons 



1 Introduction 



Since the discovery of Z)s(2317) by BaBar in April 2003, a series of new heavy 
mesons^ with open-charm and closed-charm have been observed by Belle, 
CDF, CLEO, BaBar, and BES. Among these new heavy mesons, the nar- 
row charmonium-like state X(3872) (with width < 2.3 MeV) first observed 
by Belle [5] in the exclusive decay — > K^X K^Tr'^ir^ J/ip seems to be 
the most remarkable. The evidence of X(3872) has been confirmed by three 
experiments in three decay and two production channels [6, 7, 8]. Recently, 
the decay X(3872) lJ/ip has been observed by Belle [9], which implies that 
the charge conjugation of X(3872) is positive. Further, X(3872) is unlikely 
a vector 1 state, since the bound V{e^e~)Br{X 7r"^7r~J/^) < lOeV^ at 
90% C.L. has been obtained [10], with the data collected by BES at ^/s = 4.03 
GeV. Now, the quantum numbers 0"'"+ and O""*" are ruled out based on the an- 
gular correlations in the ti^-k" J/ip system [11], also and 1 are strongly 
disfavored according to the dipion mass distribution [11]. Thus it is likely that 
X(3872) possesses J^^ = 1++ [11]. 

Theoretically, one can hardly interpret X(3872) as 2P or ID states in 
the charmonium spectrum, in view of its extremely narrow width. Thus it is 
most likely an exotic (non-gg) meson (e.g., molecule, diquark-antidiquark, and 
hybrid meson). Now the central question is whether the spectrum of QCD 
possesses a resonance around 3872 MeV with J^'-' = 1+"*". 

In this paper, we investigate the mass spectra of molecular and diquark- 
antidiquark interpolating operators whose lowest-lying states having J^'^' = 
1"*"+, in lattice QCD with exact chiral symmetry [12, 13, 14, 15, 16]. This 
study follows our recent investigation on the mass spectrum of exotic mesons 
with J^^ = 1— [17], which suggests that y(4260) [18] is in the spectrum of 
QCD, with quark content (cucu). 

For two lattice volumes 24'^ x 48 and 20'^ x 40, each of 100 gauge configura- 
tions generated with single-plaquette action at /3 = 6.1, we compute point-to- 
point quark propagators for 30 quark masses in the range 0.03 < m^a < 0.80, 
and measure the time-correlation functions of the exotic meson operators which 
can overlap with X(3872). The inverse lattice spacing is determined with 
the experimental input of while the strange quark bare mass m^a = 0.08, 
and the charm quark bare mass rrica = 0.80 are fixed such that the masses 
of the corresponding vector mesons are in good agreement with 0(1020) and 
J/'iIj{3097) respectively [19]. 

Note that we are working in the quenched approximation which in principle 
is unphysical. However, our previous results on charmed baryon masses [19], 
and also charmed meson masses and decay constants (theoretical predictions) 
[20] turn out to be in good agreement with the experimental values. This 

-'^For recent reviews of these new heavy mesons, sec, for example, Refs. [1, 2, 3, 4], and 
references therein. 
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suggests that it is plausible to use the quenched lattice QCD with exact chiral 
symmetry to investigate the mass spectra of the exotic meson operators con- 
structed in this paper, as first steps toward the unquenched calculations. The 
systematic error due to quenching can be determined only after we can repeat 
the same calculation with unquenched gauge configurations. However, the 
Monte Carlo simulation of unquenched gauge configurations for lattice QCD 
with exact chiral symmetry, on the lattices 20^ x 40 and 24^ x 48 at /3 = 6.1, 
still remains a challenge to the lattice community. Thus, in this paper, we 
proceed with the quenched approximation, assuming that the quenching error 
does not change our conclusions dramatically, in view of the good agreement 
between our previous quenched mass spectra of charmed hadrons [19, 20] and 
the experimental values. 



2 The Molecular Operators 

In this section, we construct four molecular operators with quark content 
(cqcq) such that the lowest-lying state of each operator has J^^ = l"*"^. 
Then we compute the time correlation function of each operator, and extract 
the mass of its lowest-lying state. Explicitly, these molcular operators are: 

^1 = ^ {(q7iC)(c75q) - (C7iq)(q75c)} (1) 

M2 = (q757^q)(cc) (2) 

^3 = ^ {(q757ic)(cq) + (C757iq)(qc)} (3) 

M4 = (C757ic)(qq) (4) 
The time-correlation function^ 

CMit)=Y.{M{x,t)M\0,Oy 



is measured for each gauge configuration, and its average over all gauge con- 
figurations is fitted to the usual formula 
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2ma 

to extract the mass ma of the lowest-lying state and its spectral weight 

w ' 



2ma 



^Here we have neglected the cc and qq annihilation diagrams such that C{t) does not 
overlap with any conventional meson (cc or qq) states. Also, C{t) has been averaged over Ci 
(with 7i) for i = 1,2, 3, where in each case, the "forward-propagator" Ci{t) and "backward- 
propagator" Ci{T — t) are averaged to increase the statistics. The same strategy is applied 
to all time-correlation functions in this paper. 
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Figure 1: The ratio of spectral weights of the lowest-lying state of the molecular 
operator Mi, for 20^ x 40 and 24^ x 48 lattices at /3 = 6.1. The upper-horizontal 
line R = (24/20)"^ = 1.728, is the signature of 2-particle scattering state, while 
the lower- horizontal line i? = 1.0 is the signature of a resonance. 



Theoretically, if this state is a genuine resonance, then its mass ma and spectral 
weight W should be almost constant for any lattices with the same lattice 
spacing. On the other hand, if it is a 2-particle scattering state, then its 
mass ma is sensitive to the lattice volume, and its spectral weight is inversely 
proportional to the spatial volume for lattices with the same lattice spacing. 
In the following, we shall use the ratio of the spectral weights on two spatial 
volumes 20^ and 24"^ with the same lattice spacing (/3 = 6.1) to discriminate 
whether any hadronic state under investigation is a resonance or not. 

In Fig. 1, the ratio {R = W2Q/W24) of spectral weights of the lowest-lying 
state extracted from the time-correlation function of Mi on the 20^ x 40 and 
24^ X 48 lattices is plotted versus the quark mass m^a G [0.03, 0.80]. (Here the 
quark fields q and q are always taken to be different from c and c, even in 
the limit mq — m^.) Evidently, R ~ 1.0 for the entire range of quark masses, 
which implies that there exist J^*" = resonances, with quark content 
(cscs), and (cucu) respectively. 

In Fig. 2, the mass of the lowest-lying state extracted from the molecular 
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Figure 2: The mass of the lowest-lying state of Mi versus the quark mass m^a, 
on the 24^ x 48 lattice at /3 = 6.1. The solid line is the linear fit. 



operator Mi is plotted versus m^a. In the limit — > ~ 0.00265a~^ 
(corresponding to rrij^ = 135 MeV), it gives m = 3895(27) MeV, which is in 
good agreement with the mass of X(3872). 

For rriq = rUs = 0.08a~^, the time-correlation function and effective mass 
of Ml are plotted in Fig. 3a and Fig. 3b respectively. With single exponential 
fit, it gives m[(s7jc)(c75s) - (c7js)(s75c)] = 4109(21) MeV. 

Next we turn to the molecular operators M2, M3, and M4. In each of 
these cases, the ratio of spectral weights {R = W20/W24) behaves as i? ~ 1.0 
for iriqa > 0.05, but deviates from 1.0 with large errors as — > m„. This 
seems to suggest that M2, M3 and M4 have little overlap with the resonance 
detected by Mi as nig m„. We suspect that this is due to quenching 
artifacts as rrig — > mostly coming from the scalar meson (qq), or (cq), or 
(qc), as well as the pseudovector (q757jq). However, all molecular operators 
give compatible masses as rriq This is also the case for the resonance 

at iTiq = tris- 
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Figure 3: (a) The time- correlation function C{t) of the lowest-lying state of 
Ml for ruq = ms = 0.08a"\ on the 24^ x 48 lattice at /3 = 6.1. The solid line 
is the hyperbolic-cosine fit for t G [14,21]. (b) The effective mass M^ffit) = 
ln[C{t)/C{t + 1)] of Cit) in Fig. 3a. 



3 The Diquark-Antidiquark Operator 

We construct the diquark-antidiquark operator with J^*" = 1++ as 

^4(3;) = {(q^C7iC)^„(qC75C^)^a - (qC7fc'^)^a(q^C'75c)xa} (5) 

where C is the charge conjugation operator satisfying C'-ff^C'^^ = — 7J and 
(6*75)"^ = — C75. Here the "diquark" operator (q^rQ)^.^ for any Dirac matrix 
r is defined as 

(q'^rQ)x.a = ^abc^lxab^ afsQxPc (6) 

where x, {a,b,c} and {a,P} denote the lattice site, color, and Dirac indices 
respectively, and eabc is the completely antisymmetric tensor. Thus the diquark 
(6) transforms like color anti-triplet. For F = C75, it transforms like = O"*", 
while for F = C'ji {i = 1, 2, 3), it transforms like l"*". 

In Fig. 4, the ratio {R = ^^20/^^24) of spectral weights of the lowest-lying 
state extracted from the time-correlation function of X4 on the 20'^ x 40 and 
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Figure 4: The ratio of spectral weights of the lowest-lying state of diquark- 
antidiquark operator X^, for 20^ x 40 and 24^ x 48 lattices a.t (3 = 6.1. The 
upper-horizontal line R = (24/20)^ = 1.728, is the signature of 2-particle 
scattering state, while the lower- horizontal line = 1.0 is the signature of a 
resonance. 



24^ X 48 lattices is plotted versus the quark mass m^a G [0.03, 0.8]. Evidently, 
i? ^ 1.0 for the entire range of quark masses, which implies that there exist 
J^*-^ = l"*"^ resonances, with quark content (cscs), and (cucu) respectively. 

In Fig. 5, the mass of the lowest-lying state of the diquark-antidiquark 
operator X4 is plotted versus m^a. In the limit ruq — it gives m = 
3891(17) MeV, which is in good agreement with the mass of X(3872). 

For rriq = rris = 0.08a~^, the time-correlation function and effective mass 
of the diquark-antidiquark operator are plotted in Fig. 6. The mass of the 
lowest-lying state is 4134(19) MeV. 

4 Summary and Discussions 

In this paper, we have investigated the mass spectra of several interpolating 
operators (i.e., the molecular operators Mi, M2, M3, and M4, and the diquark- 
antidiquark operator X4) with the lowest-lying J^*^ = l"*""*", in quenched lattice 
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Figure 5: The mass of the lowest-lying state of the diquark-antidiquark oper- 
ator X4 versus the quark mass m^a, on the 24^ x 48 lattice at /3 = 6.1. The 
solid line is the linear fit. 



QCD with exact chiral symmetry. Our results for Mi and X4 are summarized 
in Table 1, where in each case, the first error is statistical, and the second one is 
our estimate of combined systematic uncertainty including those coming from: 
(i) possible plateaus (fit ranges) with x^/d.o.f. < 1; (ii) the uncertainties in the 
strange quark mass and the charm quark mass; (iii) chiral extrapolation (for 
the entries containing u/d quarks); and (iv) finite size effects (by comparing 
results of two lattice sizes). Note that we cannot estimate the discretization 
error since we have been working with one lattice spacing. Even though lat- 
tice QCD with exact chiral symmetry does not have 0(a) and 0{ma) lattice 
artifacts, the 0{w?o?) effect might turn out to be not negligible for m^a = 0.8. 

Evidently, both the molecular operator Mi and the diquark-antidiquark 
operator X4 detect a l"*""*" resonance around 3890 ± 30 MeV in the limit — >■ 
which is naturally identified with X(3872). This suggests that X(3872) is 
indeed in the spectrum of QCD, with quark content (cucu), and J^^ = 

Now, in the quenched approximation, our results suggest that X(3872) has 
good overlap with the molecular operator Mi as well as the diquark-antidiquark 
operator X4. This is in contrast to the case of y(4260) in our recent study [17], 
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Figure 6: (a) The time- correlation function C{t) of the lowest-lying state of 
X4 for niq = nis = 0.08a~^, on the 24^ x 48 lattice at (3 = 6.1. The solid line 
is the hyperbolic-cosine fit for t G [14,21]. (b) The effective mass M^ffit) = 
\n[C{t)/C{t + 1)] of C{t) in Fig. 6a. 



in which y(4260) seems to have better overlap with the molecular operator 
{(q757jc)(c75q) — (c757jq)(q75c)} than any diquark-antidiquark operators. 
This seems to suggest that X(3872) is more tightly bound than y(4260). It 
would be interesting to see whether this picture persists even for unquenched 
QCD. 



Finally, for 



heavier exotic meson resonance with J^*^ 



IS 



also detected, with quark content (cscs) around 4100 ± 50 MeV, which serves 
as a prediction from lattice QCD with exact chiral symmetry. 
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Operator 


Mass (MeV) 


R/S 


71 [(u7ic)(c75u) - (C7iu)(u75c)] 
^[(S7,c)(c75s) - (C7,s)(s75c)] 


3895(27) (35) 
4109(21)(32) 


R 
R 


V2 ' 

1 

V2 


^(u^C7^c)(uC75C^) - (u^C75c)(aC7rc^)} 
{(s^C7,c)(sC75C^) - (s^C75c)(sC7fc^)} 


3891(17)(21) 
4134(19)(25) 


R 
R 



Table 1: Mass spectra of the molecular operator Mi and the diquark- 
antidiquark operator X4 with J^'^ = 1++. The last column R/S denotes 
resonance (R) or scattering (S) state. 
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